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Abstract 

Beauty and charm e~^e~ factories running at resonance thresholds have unique capabihties for studies 



o 

(N 

^ . of the production of hght Dark Matter particles in the decays of Bq[D) meson pairs. We provide a 
comprehensive study of light Dark Matter production in heavy meson decays with missing energy ^ in 

(N 

J> ' the final state, such as Bq{D^) — )• ^ and Bq{D^) — )• 7^. We argue that such transitions can be studied 

C^^ . at the current flavor factories (and future super-flavor factories) by tagging the missing-energy decays 
CN ■ 



with B{D ) decays "on the other side." 
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I. INTRODUCTION 



The presence of cold Dark Matter (DM) in our universe provides the most natural explanation 



for several observational puzzles, from the original measurement of the rotational curves [l| of 



fl 



galaxies to the observation of background objects in the Bullet Cluster [2| and spectrum features of 
the cosmic microwave background (CMB) fluctuations. In the conventional picture, DM accounts 
for the majority of mass in our Universe. However, the nature of DM is still very much a mystery, 
which could intimately connect astronomical observations with predictions of various elementary 
particle theories. Many such theories, with the notable exception of the Standard Model (SM), 
predict one or more stable, electrically-neutral particles in their spectrum |3|. These particles 
could form all or part of the non-baryonic Dark Matter in the Universe. 

Different models provide different assignments for DM particles' spin and various windows for 
their masses and couplings to luminous matter. In the most popular models DM is a weakly 
interacting particle particle with mass set around the electroweak energy scale. This follows from 



a 



the experimental measurements of the relic abundance VtoMh"^ ~ 0.12 by WMAP collaboration J4| 

VtoMh'^ ~ {cTannVrel) ^ OC -— ~ 0.12, (1) 

9 

where M and g are the mass and the interaction strength associated with DM annihilation re- 
spectively. As one can see, a weakly-interacting massive particle (WIMP) with electroweak-scale 
mass naturally gives the result of Eq. ([T]). This, coupled with an observation that very light DM 
particles might overdose the Universe (known as the Lee- Weinberg limit (5|) , seems to exclude 
the possibility of the light-mass solution for DM, setting Mdm > 2 — 6 GeV. 

A detailed look at this argument reveals that those constraints could be easily avoided, so even 



BE 



MeV-scale particles can be good DM candidates. For instance, DM could be non-fermionic |6|, |7|, 
in which case the usual suppression of the DM annihilation cross-section used in setting the Lee- 
Weinberg limit does not hold. In addition, low energy resonances could enhance the cross-section 
without the need for a large coupling constant. Other solutions, which also provide low- mass 



candidates for DM particles, are also possible |8l-ll0|. 



There are many experiments designed to search for both direct interactions of DM with the 



detector and indirect evidence of DM annihilations in our or other galaxies by looking for the 



ow- 



ul 



products such as gamma-rays, positrons and antiprotons. Those can in principle probe 
mass DM. However, direct searches, performed by experiments such as DAMA and CDMS 
rely on the measurement of the kinematic recoil of the nuclei in DM interactions. For cold DM 
particles, such measurements lose sensitivity with the decreasing mass of the WIMP as recoil 



energy becomes smaller [12|]. Indirect experiments, such as HESS [13|], are specifically tuned to 
see large energy secondaries, only possible for weak-scale WIMPs. The backgrounds for positron 
and antiproton searches by HEAT and/or PAMELA experiments [ij] could be prohibitively large 
at small energies. 

It is well-known that the existing e~^e~ flavor factories and future super-flavor factories could 
provide the perfect opportunity to search for rare processes, especially the ones that require high 
purity of the final states. In particular, probes of rare B-decays, such as 5 — )■ K^*^uV^ are only 
possible at those machines. These colliders, where Bg{D) and Bq{D) are produced in charge and 
CP-correlated states, have an opportunity to tag the decaying heavy meson "on the other side," 
which provides the charge or CP-identification of the decaying "signal" B or D meson. In fact. 



many CP- violating parameters at B-factories have been measured using this method [15|. It is then 
possible to perform a similar tag on the meson decaying to a pair of light DM particles or a pair of 
DM particles and a photon. The latter process might become important for some DM models as it 
eliminates helicity suppression of the final state^. Moreover, compared to B ^ K -\- ^ transitions, 
where ^ is missing energy, a massless photon could provide better experimental opportunities for 
tagging without reducing the probed parameter space of the DM masses. Finally, searches for light 
DM in heavy meson decays could be more sensitive than direct detection and other experiments, 
as DM couplings to heavy quarks could be enhanced, as for example happens in Higgs portal 



models |l6 |. 



In this paper we compute branching ratios for the heavy meson states decaying into XsXs ^^^ 
XsXsl- Here Xs is a DM particle of spin s, which appears as missing energy in a detector. The 
DM anti-particle Xs ^^^J o^ ^^^ly not coincide with Xs- We shall first consider model-independent 



^ This is similar to the situation in Icptonic decays of 5-niesons, where the branching ratios B{B — > ^v^) 
B{B -^ fiiy) and B{B -^ eu-f) > B{B -^ ev). 



interactions of DM particles of spin-0, spin-1/2, and spin-1 with quarks. In each case we write 
the most general effective Hamiltonian coupling DM particles to flavor- changing 6 — )■ g (where 
q = s{d)) or c — )■ M current and compute B{D) — )• XsXsil) decay rates. We then consider popular 
models, already available in the literature, that can generate those processes. 

II. FORMALISM AND THE STANDARD MODEL BACKGROUND 

The computation of decay rates for two-body processes Bg{D) — )• XsXs is a straightforward 
task which only requires the knowledge of appropriate B — t- vacuum matrix elements. We use 
conventional parameterization for those. 



(0| bf^q \B,)=0, 



(0| bq \B,) = 0, 



(0| 67^75? \B,)=zfsP>', (0| 675? \B,) 



■I 

mi, + niq 



(2) 



where P^ is the 4-momentum of heavy meson Bq. Similar formulas can be obtained for Z^-meson. 
In what follows we shall provide relevant derivations for Bq mesons only, but report results for 
both Bq and D'^-meson decays. 

Before computing the relevant DM production rates, let us study the Standard Model back- 
ground for the decays with missing energy realized in transitions to uV states. The Standard 
Model effective Hamiltonian for Bq{D) — )■ vT'{'-)) reads 
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AGt 
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y/2 27r sin^ 9 
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E Y.^kX\xu){rQq){v^LiA)^ 



(3) 



=e,/i,r k 



where Jq = qi'^^bi for beauty, and Jqq = UlI^Cl for charm transitions, and we consider Dirac 
neutrinos. The functions AfcX'(xfc) are relevant combinations of the Cabbibo-Kobayashi-Maskawa 
(CKM) factors and Inami-Lim functions. For 6 — )■ g transitions these functions are overwhelmingly 
dominated by the top-quark contribution. 



Xt 



Y,XkX\xk) = v;qVtbX{xt), with x{xt) = ^ 



' xt + 2 3{xt-2) 

xt-1 (xt - ly 



Inxt 



(4) 



and Xt = ml/M'^. Perturbative QCD corrections can be taken into account by the replace- 
ment 17 1 

r rv 1 r rv / _ 9?i \ " 

(5) 
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where Xi{xt) can be found in Ref. [l7|. They change our estimate by at most 10%, and therefore 
be neglected in our analysis. For c — ?■ n transitions we keep the contributions from both internal 
h and s- quarks, so 



Y.\uX\xk) = v:yusX' 



dj a 



V:,VubX\x,), with X\x,) = D{x,,yi)l2 



(6) 



where D{xq, yi) is the Inami-Lim function 18|] for yi = mf/m 



^Wy 



D{xq,yi)=- 
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yi \yi 



yi-^\\ ^1 



Xn 



yi -Xq 



Xq-A 
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Xa 



yi-'^J Xq 
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(7) 



Given this, one can easily estimate branching ratios for Bq{D) — )• uV decays. One can immediately 
notice that the left-handed structure of the Hamiltonian should result in helicity suppression of 
those transitions. Assuming for neutrino masses that rriy ~ J2ii^u^ < 0.62 eV 19|], where m,^. is 
the mass of one of the neutrinos, we obtain for the branching ratio 



B{Bs -^ vV) 



lOvr^ sin" 9^^ Bs 



2„2 



VtbV:,\'X{xtYx 



3.07 X 10 



-24 



(8) 



where x^ = my/MB^ and Tb^ = ^b^ = 1/tb is the total width of the Bg meson. With tb = 1.548 ps 
we obtain B{Bd — )■ uu) = 1.24 x 10~^^. A similar calculation yields S(-D° — )■ uV) = 1.1 x 10~^°. 
Clearly such tiny rates imply that decays of heavy mesons into neutrino-antineutrino final states 
in the Standard Model can be safely neglected as sources of background in the searches for DM 
in i?q (D)-decays. This is one of the main differences between this study and studies of DM 
production in i? — > K^*^ + ^ transitions |6|. 

Helicity suppression in the final state can be overcome by adding a third particle, such as a 
photon, to the final state. The calculation of B(D) — )■ i/TJ'j has been done before [20|, so here we 



simply present an update. The branching ratio for B{D) — )■ uw-f in principle depends on several 
form-factors, 

{j{k)\h,q\B,{k + q))=e e^.^^e* V^;'^ ^^ ^^ ^ 



M, 



Ba 



{-f{k)\b-f^-f5q\Bg{k + q))=-ie e* (kg) - (e*g) k 
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M 



Ba 
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{-f{k)\ha^,q\B^{k + q)) 



Mi 



"C/lvAcr 



* „\„\l 



Ge*^k'' + He*^q'' + N{e*q)q^k'' 



(10) 



Matrix element {'~^{k)\bafj,y'~^^q\Bq{k + q)) can be obtained using identity cr^i^ = — |e^jya/30"Q,^75 [21| 



^=4^1, 

A{qk) 



N 



H 



q' 



{fi + 9i), fiigi] 



'ifi+9i] 

foigo) 



gV/i' 
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where /o, go, fif, fig are known from QCD light-cone sum rules. Similar formulas hold for i5-decays. 
It is important to note that only one out of two form- factors is independent. Indeed, as it was 



shown in 
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fBaM^ 



2E^ ^' 



(12) 



where R ^ ~ Mn - m,,, and Fb„ = -QgR, 



Ql, MBqQb~mbiQb + Qq) 



Q-^^q 



mt 



m6(A/s -mt) 



• Qq = Qb = +1/3 are the 



electrical charges of q and 6-quarks. Similar form factor can be obtained for the D-meson after a 
suitable redefinition of quark masses and charges. One-loop QCD corrections to the Eq. (TT2|) can 
also be computed [24 1. 

The amplitude for Bq{D) — ?■ uV^ transition could be written as 



A{Bg -^ i^V-f) 



2eCf^{xt) 



M, 



Ba 



e;..p.e*^g^r f^{q') + z fe (kq) - {e*q) kj f^{q')] Vl^ul, (13) 



where Cf'^'{xt) = G F(yVtbXQ{xt) / {2\^tt sin^ 9w) and e is the electric charge. This results in the 



photon spectrum and a branching ratio, 



dV 4/1 GW 
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(14) 



BiB,^u,^^)=—flGWMl,jV,,V;,Xo{x,)fij^-^^^j , (15) 



where we set x^, = 0. Numerically, B{Bs — ?■ 1^2/7) = r(i?s — ?• z^i^7)/rBs = 3.68 x 10^^. Similar 
results for Bd and D^ mesons are B{Bd — )■ uu^y) = 1.96 x 10~^ and B{D'^ — )■ 2/1^7) = 3.96 x 10~^^ 
respectively. 

It is important to notice that the approach to rare radiative transitions described above works 
extremely well for SM neutrinos in the final state since E^ ^ ^qcd over most of the available 
phase space. It might not be the case for the DM production. In particular, for rrioM ^ 2 GeV, 
the photon energy is quite small and corrections to Eq. flT2|) could become significant. Therefore, 
our results obtained by using the formalism above should be corrected, for instance, using heavy 
meson chiral techniques. 

Currently the only experimental constraints on Bq{D^) — t- ^ and Bq{D^) -^ 7^ transitions are 



available from Bd decays 25 1 



B{Bd ^ ^) < 2.2 X 10-^ 

i3(5rf ^ ^ + 7) < 4.7 X 10-^ (16) 

One can see that while the branching ratios for the decays into z/z77 final states are orders of 
magnitude larger than the corresponding decays into uu final states, they are still way beyond 
experimental sensitivities of currently operating detectors. Thus, we conclude that SM provides 
no irreducible background to studies of light DM in such decays. 



III. SCALAR DARK MATTER PRODUCTION 

A. Generic effective Hamiltonian and B — t- XoXoIt) decays 

Let us consider the generic case of a complex neutral scalar field xo describing the DM and 
limit our discussion to effective operators of dimensions no more than six. In this case, a generic 
effective Hamiltonian has a very simple form, 

^5^/ = 2E^O., (17) 

i 

where A is the scale associated with the particle (s) mediating interactions between the SM and 
DM fields, and C^ are the Wilson coefficients. The effective operators are 

Oi=mb(bRqL)ix*oXo), 

02=mb(bLqR)iXoXo), (18) 

03={bLl''qL)iXod^,Xo), 

04={bRYqR){x*odf.Xo), 

•O- -)> <- 

where d= [d — 9)/2. For relevant D-meson decays one should substitute mf, — )• m^ and b ^ q 
currents with c — )■ m currents. Operators 03,4 disappear for DM in the form of real scalar fields. 
We note that while the generic form of Eq. f fTSj) implies that the mediator of interaction between 
DM and the SM fields is assumed to be heavy. Ma > rriB (d), it is easy to account for the light 
mediator by substituting C^/Is?' — )■ C^ /{M"^ rj-,^ — M^). Clearly, a resonance enhancement of 
B[D) — 7- XoXo I'ate is possible if for some reason the mediator's mass happens to be close to 
^Bq{D)- If observed, this resonance enhancement would be seen as anomalously large Wilson 
coefficients of the effective Hamiltonian of Eq. flTSj) . 

Let us first compute the B{D) — )• XoXo transition rate. It follows from Eq. (TT8|) that the decay 



branching ratio is 



BIB :rr) K' " ^H ( fB,Mi,,n, x' 



4x1 



(19) 



where x^ = m^/ Mb is a rescaled DM mass. Clearly, this rate is not helicity- suppressed, so it 
could be quite a sensitive tool to determine DM properties at e^e~ flavor factories. The result for a 
corresponding D-decay can be obtained via trivial substitution of quark masses, widths and decay 
constants. Computing the decay rate for various values of Dark Matter masses and comparing it 
with the experimental results for 5^ missing energy decays 
constraints on coupling constants: 



25| from Eq. (TT6|) we get the following 
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<2.03 X 10"^^ GeV"^ for m^ = 



<2.07 X 10"^^ GeV"^ for m^ = 0.1 x Mb, 



<2.22 X 10"^^ GeV^^ for m^ = 0.2 x Mb, 



<2.54 X 10"^^ GeV"^ for m^ = 0.3 x Mb, 



<3.39 X 10-^^ GeV"^ for m^ = 0.4 x M 



B, 



(20) 
(21) 
(22) 
(23) 
(24) 



It is worth pointing out that constraints obtained here are much stricter than those in 26 1. 

Applying the formalism described above, distribution of the photon energy and decay width of 
radiative decay Bg{D) — )• XoXo7 can be computed. 



dT.^ , , flaCi'^d:^ fFBX2MlE,{MB,{l-ixl)-2E^)y^ 



dE, 



B{B, ^ xlXol)-- 



3A4 X^Ti ) 

flaCi'^d'^Ml (Fb^ 
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Ba 
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We observe that Eqs. (125 p and (126 p do not depend on C| 2- This can be most easily seen from the 
fact that Bq{D) — )■ 7 form factors of scalar and pseudoscalar currents are zero, as follows from 
Eq. iQ. Computing decay rates for various values of Dark Matter mass we are able to restrict DM 
properties based on experimental constraints on B^ decays with missing energy given in Eq. (TT6|) : 

^^ < 1.55 X 10-^2 Q^y-i for ^ = 

A2 A2 - 

-^-^ < 1.86 X 10"^^ GeV-^ for m = 0.1 x Mb, 
A2 A2 - '"'' 

^^ < 3.20 X 10-^2 Q^y-A for m = 0.2 X Mb, (27) 

A2 A2 - '''' ^ ^ 



A2 A2 



< 9.06 X 10"'^ GeV-^ for m = 0.3 x Ms, 



cr^ a 



'^ ^ < 7.44 X 10'^^ GeV-^ for m = 0.4 x M, 



Srf 



Note that Eqs. ( 125|) and ( 126|) depend on C3 and C4, while Eq. (TT9|) only on Ci and C2. Since the 
models with self-conjugated DM scalar fields only contain operators Oi and O2, Bq[D) — )■ XoXo(7) 
transitions could be used to test the structure of the scalar DM sector. 

B. Production rates in particular models with scalar DM 

In this section we apply the techniques described above for the most general effective Hamil- 
tonian for DM particles interacting with the SM fields to particular model implementations of 
scalar DM, already available in the literature. The list of models considered below is by no means 
exhaustive. 
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1. Minimal and next-to-minimal Scalar Dark Matter models 

The simplest possible model for scalar DM involves a real scalar field xg ^ S coupled to the SM 
particles through the exchange of Higgs boson [y, \2j\ (see also 28|). This is also a very constrained 



model, where the only two new parameters are the mass parameter ttiq of the scalar DM particle 
S and the Higgs-scalar coupling A. Nevertheless, it is possible to have light DM in this model 
even though it might require some degree of fine-tuning. The SM Lagrangian is modified by 

4 2 

=^S' + ^K + A^vf)^' + ^VEwS'h + ^S^h' (28) 

where H is the Standard Model Higgs doublet, vew = 246 GeV is the Higgs vacuum expectation 
value and h is the corresponding physical Higgs boson. We require S to satisfy S -> —S to 
make it a good Dark Matter candidate. The scalar DM particle can be made light by requiring 
cancellations between the terms defining its mass, m? = rn^ + \v%^r. 

The transition B — )■ SS occurs in the minimal model as a one-loop process, and since mediating 
Higgs boson is much heavier than other particles involved in the process, it can be integrated out. 
The resulting effective Hamiltonian reads 



,{s) 



?>\glVtsV;^Xtmb 



which implies that c[%^ = 0, C^'^ = 3\glVtsVtlxt/128n^, and A = Mh- Thus, from Eq. ((19]), the 
branching ratio for the B — )■ SS decay in this model is 



B{B^ -^ SS) 



1287r2 



4x1 f \'\f IbM. 



IGttMbTb, \MfjJ \mb 



(30) 



where Xs = ms/niB ■ Note that this rate depends not only on the mass of S but also on the 
parameter k = X^/Mfj. This parameter also drives the calculation of the relic density of S [27|, 



-^ — X lim — ^ , 



(TannVrel = TTg ^ 1™ ^' (31) 
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FIG. 1: B{Bd — >■ SS) as a function of x = itls/Mb^- Values of A and M^ were fixed at 1 and 120 GeV 
respectively 

where Th*x is the rate for the decay h* ^ X for a virtual Higgs with Mh ~ '^f^s- We can, 
therefore, fix k from the rehc density calculation. This gives for the branching ratios of Bg and 
D-decays, 



\2 

B{B,, -^ SS)^{a.5 X 10^ GeV^) x —^Jl- 4x% 

A2 



B{Bd -^ SS)^( 1.3 X 10^ GeV^ x 4^\/l - 4x1 



Mh 



B(D° -^ SS)^(2.9 X 10-^ GeV^) x 



M'h 



l-4x? 



(32) 
(33) 
(34) 



We require the branching ratios to be smaller than the current experimental upper bound 25| for 
the missing energy decay given in Eq. ( 1T6|) . With this we are able to put the following restriction 
onto the parameters of this model: 



A 



M, 



l-4x| < 1.68 X 10"^. 



(35) 



We present the resulting branching ratios as a function of m^^-, in Fig. [H Comparing the above 
branching ratio with the available experimental data we can put constraints on the parameters of 
this model, which we present in Fig. [2l For the particular values of Dark Matter particles mass 
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FIG. 2: (a) allowed values of the DM-Higgs coupling A as a function of x = ras/Ms^ (below the curves) 
for the Higgs masses of 110 GeV (red), 120 GeV (green), and 150 GeV (blue), (b) Allowed values of 
the Higgs mass in GeV (above the curves) for A = 0.1 (red), 1 (green), and 5 (blue) as a function of 



we get 



A 



A 



Ml 
A 



Ml 
A 



A 



M2 



H 



<8.2 X 10"^ GeV"^ for m^ = 
<8.3 X 10^'^ GeV-2 for m^ = 0.1 x Mb^ 
<8.6 X 10~'^ GeV-2 for ms = 0.2 x Ms, 
<9.2 X 10-"^ GeV-2 for ms = 0.3 x Mb^ 
<1.1 X 10"^ GeV-2 for rug = 0.1 x Mb 



(36) 



The minimal scalar model described above can be made less restricted if we introduce another 
mediator for DM-SM interactions, which should somewhat alleviate the fine-tuning present in the 
minimal model 27|]. This can be done in a variety of ways. The simplest one is to introduce 
another Higgs-like field U, 



2 2 

= ^S^ + ^M^ + fiuS"^ + riVEwuh + . . . , 



(37) 
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where we only display mass and relevant intreaction terms; ellipses stands for other terms in the 
Lagrangian that are irrelevant for this discussion. 

Here u denotes the excitation around vev of U, and fi and i] are parameters with values of 
the order of electroweak scale. As far as the studies of DM production in heavy flavor decays 
are concerned, extended models of this class are equivalent to the minimal model after suitable 



redefinition of parameters 27| . Performing such redefinitions, we obtain 



2, ,2 



BiB, ^ SS)^{2.1 X 10-^) X 1^^1-Axl 

2 2 

B{B, -^ SS)^{6.3 X 10-6) X !L^^l-Axl (38) 



M^ 



n2,,2 



BiD' ^ SS)^{1.38 X 10-1^) X ^^1-Axl 

where Mu is the mass of the Higgs-like field U of Eq. ( 137|) . In the results above, the mass of the 
Higgs boson was fixed at Mh = 120 GeV. Since the S'-field is a real scalar field in both the minimal 
and the extended models, these models do not give rise to the radiative decay Bg — t- SS^. 

2. Dark Matter with two Higgs doublets (2HDM) 

In this subsection we consider a singlet scalar WIMP S that interacts with two Higgs doublets, 

-C=^S' + \^S^\H',\' + \H^\') + X2SMh:\^ + IH^]^) + hS^HjH^ - H'.HI). (39) 

We shall assume that Ai ^ A2, as the opposite limit gives results that are not different from 
the minimal scalar model considered above. The contribution of A3 is suppressed because of the 
cancelation of two diagrams, as explained in |6| . 

Calculating the effective Hamiltonian results in the following expressions of the Wilson coeffi- 
cients, 

C, -C, i2M\l-a,Y andA-M^, (40) 

where a^ = {mq/MuY. As in the previous subsection, no decay into the dark matter with photon 
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possible within the framework of this model. However, decay into a pair of dark matter particles 
is possible 

BiB. -, SS)^{0.73 X 10^ GeV^) x A^y^^^ ( " M^l - a^P ' )' ' 

^atlogat -at + 1^ 



B{B, ^ 55)^(2.1 GeV^) x A^l - 4x| [ ^'^1(1-2)'^ ' ' ^^^^ 



S(D° ^ SS)^{5.0 X 10^ GeV^) x A^^T^^ ( ^ ^^.^4^^ 



2 

* Og log Og - Og + 1 



ag)2 



Eqs. ( 14T|) can be used for constraining parameters of this model in Bq — ?> SS transitions. 

IV. FERMIONIC DARK MATTER PRODUCTION 

A. Generic effective Hamiltonian and Bm^-^ — )• Xi/2Xi/2(7) decays 

Let us now consider a generic case of fermionic Dark Matter production. It is possible that 
the DM particles have half-integral spin; so many New Physics models, including Minimal Su- 
persymmetric Standard Model (MSSM), have fermionic DM candidates. Most of those models, 
however, naturally assign rather large masses to their DM candidates. Nevertheless, either after 
some fine-tuning of the relevant parameters or after introducing a light DM-SM mediator, rel- 
atively light DM particles are still possible. Let us consider their production in the decays of 
heavy mesons. Once again, limiting ourselves to the operators of dimension of no more than six, 
a relevant effective Hamiltonian reads 

^f}/ = ^E^F^Q^> (42) 

i 

where Cj's are relevant Wilson coefficient and A represents the mass scale relevant for DM-quark 
interactions (e.g. mediator mass). In general, there are twelve possible effective operators, 

Ql=(bLl^,SL){Xl/2L^^Xl/2L), Q2 = (5l7/.Sl)Xi/2^7''Xi/2r), 
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Q3=(bR'y^SR){Xi/2^YXl/2L)^ Qi= ibR-ff,SR){Xl/2j{f''Xl/2ji), 

Q5 = (bLSR)(Xl/2LXl/2R), Qe = {bLSR){Xl/2RXl/2L), (43) 

<57=(5rSl)(Xi/2^Xi/2r), Qs = (PrSl){Xi/2rXi/2j), 

Q9 = {bL(r^,uSR){Xl/2L^'"'Xl/2R), QlO = (^Lf^M^Sfl) (Xl/2^C^^''Xl/2i), 

Qn = (bR(^,^,ySL){Xl/2L^'"'Xl/2R), Ql2 = (Pr(^iiuSl){Xi/2r(^^'' Xl/2i)^ 

where the Dark Matter fermion xi/2 can be either of Dirac or Majorana type. The latter choice 
leads to some simplification of the basis. All needed matrix elements have been given in Eq. ([2]). 
Note that the matrix elements of the tensor operators vanish, 



{Q\ha^''PL,Rq\B,) = 0. 



(44) 



For relevant D-meson decays one should substitute mi, — )■ ttIc and h ^ q currents with c — )■ u 
currents. Using the Hamiltonian of Eq. (143|) we get for the branching ration of Bq — )■ X1/2X1/2, 



B{Bq -^ X1/2X1/2) 



p Ml , 



CpsyCi 



57'-^68" 



4M1^4 
rrib + ruq) 



(C57 + Cgg) 



■2Ci_ 



! \ 

mi, + m. 



+ 2(Ci3 + C24)'x2 



^1(24-; 
(mfe + mqf 



(45) 



if) 
where we employed short-hand notations for the combinations of Wilson coefficients Cij = CI — 

( f) ~ 

Cj , and Ci-s = C13C57 + C24C57 + CisCes + C2aCqs,. Due to its larger mass chirality suppression 

for the GeV-scale Dark Matter is not as severe as for neutrinos, even for purely left-handed 

interactions. The obtained result leads to model-independent constraints on the Wilson coefficients 

of Eq. (I42p . which is based on experimental data for missing energy decays of B^ meson (see, e.g. 

Eq. ( TT6|) ). They can be found in Table [H The results presented there can be used to constrain 

parameters of particular models of fermionic Dark Matter considered below. 

The technique which we use for the computation of T{Bq{D) — > X1/2X1/27) is very similar to 

the one used for the radiative decay of heavy meson into scalar DM particles discussed above. The 

hadronic part of the matrix element remains the same, we only modify the part that describes 
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x^ 


Cl/A^ 


C2/A^ 


C3/A^ 


C4/A^ 


C,/A', 


Ce/A', 


Cr/A', 


Cs/A^ 




GeV-2 


GeV-2 


GeV-2 


GeV-2 


GeV-2 


GeV-2 


GeV-2 


GeV-2 





- 


- 


- 


- 


2.3 X 10-8 


2.3 X 10-8 


2.3 X 10-8 


2.3 X 10-8 


0.1 


1.9 X 10-^ 


1.9 X 10-^ 


1.9 X 10-^ 


1.9 X 10-^ 


2.3 X 10-8 


2.3 X 10-8 


2.3 X 10-8 


2.3 X 10-8 


0.2 


9.7 X 10-8 


9.7 X 10-8 


9.7 X 10-8 


9.7 X 10-8 


2.5 X 10-8 


2.5 X 10-8 


2.5 X 10-8 


2.5 X 10-8 


0.3 


6.9 X 10-8 


6.9 X 10-8 


6.9 X 10-8 


6.9 X 10-8 


2.8 X 10-8 


2.8 X 10-8 


2.8 X 10-8 


2.8 X 10-8 


0.4 


6.0 X 10-8 


6.0 X 10-8 


6.0 X 10-8 


6.0 X 10-8 


3.6 X 10-8 


3.6 X 10-8 


3.6 X 10-8 


3.6 X 10-8 



TABLE I: Constraints (upper limits) on the Wilson coefficients of operators of Eq. ()43p from the Bg 
X1/2X1/2 transition. Note that operators Qg — Q12 give no contribution to this decay. 

Dark Matter. These lead to 



dr dVi^s rfr9„i2 



dE^ dE^ 



dE^ 



dT^^sJiFjc^MlE, y/MB (1 - AxD - 2E, 
dE^ 247r2A2 ^Mb, - 2E^ 

X \CI + Cl + Cl + CI){Mb, - xIMb, - E^) - (3C1C2 + 'iC^C^)xlMB 

64a / 



dV. 



9-12 



^^ 



Mr ri - 4x2) - 2E, 



dE^ 3M47r2A2 \Mb, - 2E, J I Mb, - 2E^ 



X 



2 { (Cio + 9C11C10 — 3C12C10 + Cii — 3C12C11 + 3C9(Cio + Cii + Ci2))/i 



-gifiiCf, + 3Cio(Cn + C12) + CniCn + SC,^) - SC^iCo + Cu + C12)) 



+2gi{Ct, - QCoCu + Cf,)) x{M'b, 

Hfl - 9ih + 'igDiCl + Cl,){Ml^ - 2MbE,] 



(46) 



(47) 



While there are many models of light fermionic DM that employ operators Qi - Qs, we are not 
aware of the models with operators Qg - Q12. Therefore, we chose not to provide a closed analytic 
expression for Bg_i2{Bg — )■ X1/2X1/27) here due to overall bulkiness of the resulting expression. An 
interested reader can perform numerical integration of Eq. (H7|) for a particular model, if needed. 
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x^ 


Ci/A^ GeV"2 


Cs/A^, GeV-2 


Cs/A^, GeV-2 


C4/A2, GeV-2 





6.3 X 10^^ 


6.3 X 10-^ 


6.3 X 10-7 


6.3 X 10-7 


0.1 


7.0 X 10-^ 


7.0 X 10-^ 


7.0 X 10-7 


7.0 X 10-7 


0.2 


9.2 X 10~^ 


9.2 X 10-^ 


9.2 X 10-7 


9.2 X 10-7 


0.3 


1.5 X 10^6 


1.5 X 10"^ 


1.5 X 10-6 


1.5 X 10-6 


0.4 


3.4 X 10^6 


3.4 X 10-6 


3.4 X 10-6 


3.4 X 10-6 



TABLE II: Constraints (upper limits) on the Wilson coefficients of operators of Eq. (j43p from the Bq 
X1/2X1/27 transition. Note that operators Q^ — Qs give no contribution to this decay. 

Integrating Eq. ( H6|) over the photon energy analytically we obtain 



F2 /I Ml a 






(Cl + Cl + Cl + Ct) Y{x, 



9 



(C1C2 + C3C4) Z(a;J 



(4J 



where the factors Y{x^) and Z{x^) are defined as 



r(a;^) = l - 2x' + 3x^(3 - 6x' + 4x^) Jl - 4^ lo, 



Z{x^)=x;^ 



l + 2x'+8xi{l-x')Jl-Axl\og 



2x. 



1 - 4x2 ^ 



2x. 



l + ,/l-4x2^ 



llxj + 12x6, 



5X, 



(49) 



This equation can be used to place constraints on the individual Wilson coefficients of Eq. f l43|) . 
They can be found in Table HIl Both Eq. (145|) and Eq. (l48ll can now be used to constrain the 
parameters of the particular models of fermionic DM. 



B. Production rates in particular models with fermionic DM 



1. Models with hidden valleys 



It was pointed out in 30|] that there could be light particles called u-quarks interacting with 



Standard Model sector via heavy mediator Z'. In the simplest f-Model, a SU{ny) x f/(l) gauge 
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group with couplings g' and g^ is added to the Standard ModeP. The U{\) symmetry is broken by 
vacuum expectation value of the scalar field (</>), giving Z' a mass of ~ 1 — 6 TeV. The latter can 
mix with Standard Model Z via kinetic mixing kF'^'^F'. In this model the role of Dark Matter 
is played by the t^-quarks (xi/2 = v). 

The model corresponds to the following set of parameters for the decay of Bg meson (for decays 
of Bd and D^ parameters will be similar): 

'^^=^ AT ■ 2/1 ytbVt,X{x), and A = Mz' (50) 

where k is the kinetic mixing parameter, g' is a gauge coupling of the Z' and w-quarks, and Mz' 
is the mass of the heavy mediator. The rest of the Wilson coefficients Cj are zero. Thus, from 
Eq. (gSD, 

B{B, ^ vv) ^ (1.76 GeY')xl^l-Axl (j^^ (51) 

where x^ = rriy/MB ■ The corresponding results for B^ and D^ decays are 

B{B, ^ vv) ^ (4.68 X 10-^ GeV^)xl,Jl - 4x1 (j^) , (52) 

and 

g'kV 



B{D° -^ vv) ^ (2.68 X 10~^ GeV^)xl^l - 4x1 ( " 



(53) 



Mz'J 

respectively. The corresponding expression for the decay into two f-quarks and photon can be 
obtained by defining 

Gpkg'aMzMz' , .e 

^1 = o AT • 2^ — ^tbV^,X{x)-, and k = Mz'. (54) 

2(7v2 sm 6w o 



We present our results in Fig. 3(a) in order to extract the dependence on DM mass. The analytic 



The g' coupling constant introduced here is not to be confused with the SM hypercharge coupUng constant. 
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0.1 0.2 0.3 0.4 0.5 



(a) B vs m 




0.1 0.2 0.3 0.4 0.5 



(b) 



FIG. 3: (a) B{Bd — >• vv) as a function of x = m^/MB^ evaluated ai g' = \^ k = 1 and Mz' = 1 TeV; (b) 
Allowed values of the M^/ mass in GeV (above the curves) for gik = 1 (black), 0.1 (red), and 10 (green) 
as a function of x = rriy/MB^- Solid lines represent the constraints from the 2-body, and the dashed ones 
- from the 3 body (radiative) decay. As one can see, the constraints on the mass of Z' are very loose. 



results for the branching ratios can be well approximated by the following formulas, 

B{B, ^ vv^) ^ (2.76 X 10-^ GeV^^^'^' - ^^"'^^ 



Mh sll - Axl 



(55) 



for the branching ratio of B^ radiative decay and 



B{Bd -^ vv-i) ^ (9.07 X 10^^ GeV' 



6 r..v2^^ X ^^^""^ 



B{D° -^ vv-f) ^ (3.68 X 10-^2 GeV^ 



gie Y{x,) 



X 



Ml J I - 4x2 



(56) 

(57) 



for Bii and D^ decays, respectively. The structure function Y{x) appearing in this equation was 
defined in Eq. 



2. Right-handed massive neutrinos as a Ferniionic Dark Matter 



Massive right-handed neutrinos appear naturally in left-right symmetric models (see for exam- 
ple 311 )• "^^^ see-saw mechanism is used to get light left-handed neutrinos and massive right- 
handed ones. The coupling of the massive neutrino to the SM fields in this case is mediated by a 
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right-handed gauge boson with mass in the TeV range. In this section xi/2 = J^r- 

^eff = ^ n^,:^2n E ^kX{Xk) (j^,) (VR-f^UR) , (58) 

V2 zTT sm aw k 

where Jq = QR^^hR for beauty and Jq^ = ur'-j^cr for charm transitions. The functions AfcX(xfc) 
are the combinations of the Cabbibo-Kobayashi-Maskawa (CKM) factors and Inami-Lim functions. 
Gp is defined similarly to the usual Fermi constant, 

which implies that 

C4 = ^ ^^-— . (60) 

8 27r sin^ ^,^ ^ ' 

Following the procedure described above, we obtain the following results for decay branching 
ratios, 

X 3.6 X 10^ GeV^ 2 I TT / x 

B{B, ^ urUr)^ -3 xl^l - Axl , (61) 

B{B, ^ z/H.^«7)^ °'^L?'^' X r(x.), (62) 

lOH^eV^ 

M,1 



B{B, ^ z/,i^fi)^ ^\,7' x^y'l - Axl (63) 

X 1.9 X 10-2 GeV^ ,,, , 

i3(E, ^ z/flz/fi7)~ ^Ti X ^(^<^). (64) 

^/^n X 5.6 X 10-5 GeV^ 2 / —7 , , 

B{D' ^ z/^i^^)^ -p x^yi - 4x2, (65) 

«mo - A 7.6 X 10-9 GeV^ ^, , ,^^, 

i3(Z^" ^ urUr-i)^ —^ X F(x^), (66) 

where Y{x) is defined in Eq. (H9l) . These results are also presented in Fig. HI 
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1. xlO"" 
1. xlO"'^ 
1. xlO"" 
1. xlO"" 
1. xlO"^' 
1. xlO"-' 



n 



I 

I 0.01 

o!l (K2 otl (kI ots " 



0.2 0.3 0.4 0.5 



(a) 



(b) 



FIG. 4: (a) B^Bd — )• vrvr) as a function of x = niup^/MB^ evaluated at Mwr = 1 TeV , (b) Allowed 
values of the Mwr mass in GeV (above the curves) as a function of a; = rriuji/MB^- Solid lines represent 
the constraints from the 2-body, and the dashed ones - from the 3 body (radiative) decay. As one can 
see, the constraints on the mass of Wr are very loose. 

C. Major ana fermions 



Majorana particles Xi/2 = X often appear in many models of physics beyond the Standard 
Model. For generic studies of decays of heavy mesons to Majorana DM particles we can also use 
Lagrangian of Eq. ( l43l) . The resulting formulas, however, will be simplified due to the known 
properties of Majorana fermions J32|, 



Taking into account the conditions of Eq. (l67l) . we can obtain the branching ratio for Bq — t- xX 
decay, 



B{B, ^ XX)-- 



PbM 



/I _ 4-7^2 

167rrB,(m, + m,)A2V x 



'-^57 + '-^68 " 2a; ((^57 — Cgs)" 



The photon energy distribution in Bq — )■ XXI decay reads 



(67) 



dV flFlaMlE,^MB,il-Axl)-2E,^^^^^ 



dE^ 487r2A2 



Mb, - 2E^ 



X {Ct, + Ci,){MB,{l + 2x') + E,), (68) 



22 



which can be integrated over to obtain the branching fraction 

B{B, ^ XX7)= ;^3^;,^,M Ci2 + CD X (69) 

364 i°g 1^ ^2- ^ ^ + (4 + ^K + K)^^^^l 



E 



As an example, we consider a reahzation of the fermionic dark matter scenario proposed in (6|. 
In this model the Majorana fermion coupled to a higgs-higgsino pair is considered. It must be 
noted that by "higgsino" we mean a fermionic field with the same quantum numbers as a Higgs 
field. We, however, do not place any supersymmetric requirements on the coupling constants. 
With that, 

M - = ~ 

- Cf=—i^i^ + liHdH^ + Xd^HdHd + Xu^H^H^, (70) 

where M '^ fi, XuVu- The Dark Matter candidate is the lightest mass eigenstate, which we define 

as 

\ 2 2 

;^=-^cos^ + i^rfSin^, sin^^- "^" 



mi=M 1 



Xlvl ^ 



u u 



We are thus led to the following effective Lagrangian, 

where at = mf/Ml and tan/3 = Vu/vd- Matching this Lagrangian to Eqs. (H2| |43|1 . we observe 
that C5 = Cq, and the remaining coefficients Cj = 0. In addition, 

VtsVtltanb ( XaXuVufJ \ nibml In at 
Taking into account Eq. ( l69l) we conclude that no decay into xx7 is possible in this particular 
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1. xlO"' 

5.X10-" ■ 

l.xlO-" ■ 
5. xlO"" ■ 



1. xlO"" 
5. xlO"'^ 



1. xlO"'^ 



0.1 0.2 0.3 0.4 0.5 



0.1 0.2 



0.4 0.5 



(a) 



(b) 



FIG. 5: (a) B{Ba — )■ xx) a-s a function of rr = itl^/Mb^- The following numerical values were used: 
K = {XdXuVulJ')/i>Hi'vl + Ai^) = 1, tan^ = 10, Mh = 102 GeV (b) Allowed values of the k (above the 
curves) for the values of of tan/3 = 1 (red), 10 (green), 100 (blue), and 1000 (purple) while mass of Higgs 
boson was fixed at M^ = 120 GeV as a function of x = m^/MB^- 



model. However, a simpler decay into xX is possible. 



B{Bs -> XX)~l-47 X 10-^^/1 - 44 



B{Bd -^ XX)~4.16 X 10~'\/1 - Ax 



log^(Qt) / tan(/3)^„ArfA„/i ' 

(i-«*)n KK + ^^') , 

- log^(at) / tan(/3)t;„AdA„/i ' 



S(DO -> XX)-1.81 X 10--yrTi4 Aan(/3). A.A ,. 



)2 



(73) 
(74) 

(75) 



where a^ = (mg/M//)^ and x^ = m^jMB^- These results can be used to constrain the parameters 
of this model. 

V. VECTOR DARK MATTER PRODUCTION. GENERIC EFFECTIVE HAMILTO- 

NIAN AND 5q(D°) -^ XiXi DECAYS 

Vector DM is a quite popular concept in non-supersymmetric solutions of the hierarchy problem. 
In particular, it can be encountered in models with Universal Extra Dimensions (UED), little 
Higgs models with T-parity, and some variations of Randall-Sundrum models. All of the proposed 
models that the authors are aware of involve weak-scale DM particles. This however, does not 
preclude the existence of the low mass vector DM. 

Let us consider a generic case of a vector field Xi describing Dark Matter. This DM particle 
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could be either a gauge boson, corresponding to some abelian or non-abelian gauge symmetry 
broken at some higher scale, or some composite state. The only assumption that we shall make is 
that Xi is odd under some Z2-type discrete symmetry, Xi ~^ ^Xi- This condition results in the 
pair-production of DM particles. 

We shall limit our discussion to the effective operators of the dimension no more than six. 
Since no gauge symmetry related to Xi is present at the scale mg, the most general effective 
Hamiltonian should be built out of the vector field Xi ^"^^ ifs field strength tensor Xi^ ■ Iii this 
case, an effective Hamiltonian has a very simple form, 

^S^/ = E^O., (76) 

i 

where A is the scale associated with the mass of the particle mediating interactions between the 
SM and DM fields, and Cl are the Wilson coefficients. The effective operators are 

Ol=mb(bLqR)Xl^,Xl, O4 = (J>Rl^qR)xTxiu^ 

02=mb(bRqL)Xl^,Xl, O5 = {bLlf^qL)Xi''xiu, (77) 

03={bLlt,qL)xTxi,y, Oe = (hR-i^qR)xTxiv, 

where xT = (V2)e^'^"'^Xio/3 ^-^d q = s,d. As before, the Hamiltonian relevant for charmed meson 
decays can be obtained by the proper substitution of 6 — )■ g current with c —> u current. 

The Bq{D) —7- XiXi transition rate can be computed using Eq. ( !77|1 . Using the form-factors 
defined in Eq. (|2]), we obtain 



flMmlJM^il-Axl 

+ (m, + m,f [^Cj, (1 - Axl) + SCf^) xj (7J 

+ 2Ci2C34(mb + niq) (l + 2xj) xj 

where Cik = Cl —Cl" and xdm = f^x/'^Bq- If is necessary to point out that Eq. (17S1) is divergent 
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•^x 


Ci/A2, GeV-2 


C2/A^ GeV-2 


C3/A^ GeV-2 


C4/A2, GeV-2 


C5/A2, GeV-2 


Ce/A^, GeV-2 











1.4 X 10-8 


1.4 X 10-8 


8.9 X 10-9 


8.9 X 10-9 


0.1 


1.2 X 10-9 


1.2 X 10-9 


1.5 X 10-8 


1.5 X 10-8 


9.1 X 10-9 


9.1 X 10-9 


0.2 


5.1 X 10-9 


5.1 X 10-9 


1.5 X 10-8 


1.5 X 10-8 


1.0 X 10-8 


1.0 X 10-8 


0.3 


1.3 X 10-^ 


1.3 X 10-*^ 


1.6 X 10-8 


1.6 X 10-8 


1.2 X 10-8 


1.2 X 10-8 


0.4 


2.9 X 10-^ 


2.9 X 10-*^ 


1.9 X 10-8 


1.9 X 10-8 


1.9 X 10-8 


1.9 X 10-8 



TABLE III: Constraints (upper limits) on the Wilson coefficients of operators of Eq. ()77p from the 
Bq — 7- XiXi transition. 



at m^ = 0, which is related to the fact that operators in Eq. ( J771) contributing to the effective 
Lagrangian are not gauge invariant. Thus, for the case of massless DM the upper Umit on the 
Wilson coefficients Ci and C2 is zero (see Table ITn|) . 

Using Eq. ( jTHj) . we can place general constraints on the Wilson coefficients of the effective 
Hamiltonian describing interactions of vector DM with quarks (see Eq. (1761) ). They are presented 
in Table Iml 

We are not aware of particular models of light DM with spin-1 particles and masses m^ < 3 GeV. 



VI. CONCLUSIONS 



We have argued that missing energy decays of the heavy mesons - B^, B^ and D^ - provide 
an important way to probe different properties of Dark Matter. Consideration of different decay 
modes - two body decays, radiative and light meson + DM decays - restricts different regions of 
the Dark Matter parameter space. Combined constraints obtained from different decay modes of 
various heavy mesons provide indispensable probe of physics beyond the Standard Model in general 
and the nature of the Dark Matter in particular. For instance, observation of Bq{D^) -t- 7^, but 
non-observation of Bq{D^) — )■ ^ transitions directly point to non-self- conjugated nature of scalar 
DM. 

We reported general constraints on the Wilson coefficients of the effective operators describing 
interactions of DM with quarks (see Tables I-III). Restrictions obtained in our paper are much 
stricter than constraints from single decay modes. Our results combined with constraints from 
astrophysical observables (for example 26|]), direct detection of Dark Matter and invisible decays 



26 



of heavy hadrons 33| could provide a full set of tools needed to test (or rule out) the models of 
light Dark Matter. 
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